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We study the interplay of spin, orbital and lattice degrees of freedom in a one-dimensional Kugel- 
Khomskii model coupled to phonons. In the vicinity of the dimer point we analyze the excitation 
spectrum, mapping the spin and orbital degrees of freedom to bond operators. In particular we 
study the renormalization of the phonon propagator due to coupling to the orbital and magnetic 
excitations. Considering both, ferro- and antiferro-orbital ordering we will show that the appear- 
ance of orbiton shake-up processes in the phonon spectrum is sensitive to the type of Jahn- Teller 
distortion. The relevance of our results for optical spectroscopy on low dimensional transition metal 
compounds will be discussed. 



I. INTRODUCTION 

Orbital ordering in three dimensions (3D) is a well es- 
tablished phenomenon in a variety of transition-metal 
compoundsi. Very recently, quasi one-dimensional 
(ID) systems with potentially similar behavior, such as 
Na2Ti2Sb20^ and NaTiSi20g^ have attracted consider- 
able interest. In these systems the coupling of orbital de- 
grees of freedom to those of the spin and probably also to 
the lattice have been suggested to lead to a rich variety of 
possible phases and non-trivial elementary excitations^. 
In principle the orbitally ordered state will allow for ele- 
mentary excitations of orbital nature, i.e. orbitons. How- 
ever, while their existence has been predicted more than 
two decades ago^ it has been only very recently that an 
indirect observation of orbitons in 3D has been claimed 
in shake-up side-bands of Raman-active phonon spectra 
in LaMnOg^. This observation underHnes the relevance 
of a simultaneous description of orbital and lattice de- 
grees of freedom. Unfortunately however, for quasi ID 
materials this is still an open issue. It is the aim of this 
work to investigate the phonon spectra of a ID model 
of mixed orbital and spin degrees of freedom by inves- 
tigating its harmonic lattice dynamics in the vicinity of 
two specifically chosen Jahn- Teller (JT) distortions. The 
outline of the paper is as follows. In section I we discuss 
the ID spin-orbital Hamiltonian and its excitations us- 
ing a bond-operator method close to the so-called dimer 
line. Sections II and III will focus on the phonon-spectra 
which form in the presence of a frozen-in, either ferro- 
or antiferro-orbital JT distortions. Conclusions will be 
presented in the final section IV. Technical details are 
deferred to appendices A through D. 



II. BOND-OPERATOR DESCRIPTION OF THE 
ID KUGEL-KHOMSKII MODEL 

In this section, and prior to discussing their effect on 
the lattice dynamics we develop a description of the or- 
bital and spin excitations. In the limit of strong Coulomb 



repulsion, orbital degeneracy in the Hubbard model can 
be described in terms of a pseudo-spin model, as was 
shown several decades ago by Kugel and Khomskii^. In 
ID and for a two-fold orbital degeneracy at 1/4-filling 
this leads to the Hamiltonian of a magnetic spin-1/2 
chain coupled to an orbital pseudo-spin- 1/2 chain by bi- 
quadratic interactions 

n 

where, apart from the spin-1/2 operator Sn, a pseudo 
spin-1/2 variable T„ with = ±1/2 corresponding to 
the orbital degrees of freedom is attached to each site 
'n\ Summation over the spin indices a = x,y,z and 
7 = x,?/,z is implied. In the remainder of this paper 
energies will be given in units of K. 

The quantum phase diagram of d has been the sub- 
ject of intense research^iSiSiifliil. For Ji, J2 > 1/4 its 
spectrum is believed to be massive^. At Ji=J2=l/4 
and A=B=Q the model has SU(4) symmetry and is 
integrablei. Along the 'dimer line' Ji = (3-|-A)/4, J2=3/4 
with —2 < A < 00 it displays a two-fold degenerate, fully 
dimerized ground state. At the 'dimer point' Ji=J2=3/4 
this ground state can be expressed as a direct matrix- 
product of alternating spin and orbital singlets with an 
energy gap of A ^ 0.375i2iil. Exactly on the dimer 
line the elementary excitations are massive solitons, i.e. 
domain walls between the two-fold degenerate ground- 
states. However, off from the dimer-line the latter readily 
bind to form generalized gapful triplet-modeA 

Prior to including the effects of a JT-distortion 
we now reformulate Q in terms of bosonic bond- 
operatorsiSiiiid. This allows for an approximate descrip- 
tion of the formation of the spin and orbital singlets in 
the ground state of the dimer phase 
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Figure 1: Dimer representation of the spin-orbital chain, n 
labels the dimer, the index 1, 2 the site on a dimer. On each 
site on the upper (lower) chain sits (pseudo-spin 
Ti,„ i.e. orbital degree of freedom) respectively. 



where s{a)n ^ and i(r)l ^ create/destroy spin(orbital) sin- 
glets and triplets and the index 1(2) labels the left (right) 
site of the dimer n in Fig. Bose operators labeled by 
roman letters refer to spin space, while greek ones refer 
to orbital space. The bond-operators have to satisfy 



s{'^)is{<j)n+t{T)lj{T)c 



(3) 



which impHes the constraint of no double-occupancy of 
the bonds by singlets or triplets. 

Inserting JSJ into (QJ an interacting Bose gas is ob- 
tained. In the limit of strong dimerization one may pro- 
ceed with the latter by invoking the Holstein-Primakoff 
(LHP) approximationiSiii. I.e. the singlets are re- 
moved using (EJ and assuming s^^-* and cr^^^ to be C- 
numbers. This is consistent with a condensation of the 
singlets in the ground state, i.e. the formation of a 
dimer state. The procedure leads to square roots of 
type s(t)(CT(t)) = [1 _ i(^)t ^^i(^)^^^]i/2 approx- 
imately, may be linearized to retain triplet interactions 
up to quartic order. Since we are interested in the effects 
of mutual interactions between the orbital and magnetic 
sector we simplify even further, by keeping those inter- 
actions which mix the orbital and magnetic triplets, but 
discarding quartic terms of purely orbital or magnetic 
nature. After some algebra we find a quadratic part of 
the Hamiltonian 



N 

Hq = --[3{Ji + J2) + J2A] 
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(4) 



with <!>„ = iti^n:ta,n,T^^yn,Tx,y^n,Tln,T,^n) and again 

a = x,y, z. n and m label dimers and i and j the compo- 
nents of $* and $. The matrix elements M^^^ are given 
in Appendix El We note, that on the quadratic level no 
mixing occurs between the spin and orbital triplets. For 



the quartic spin-orbital interactions we get 
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(5) 



h.C 



with (3 = x,y. The quadratic part Hq can be diagonal- 
ized by Fourier- and BogoHubov transformation leading 
to dispersions 



] 

t^x.k = Xl,x[X2.x + X3,xC0s{k)y' 



(6) 



(7) 



where x = a,c label new quasi-particles of purely spin 
(a) and orbital (c) nature with = Wfeajj + Vkaa,-k 

and f, = 9kc\ k + ^kCa,-k- Explicit expressions for the 
dispersions u>x.k, the Bogoliubov coefficients, and Xi.x are 
listed in Appendix El After Fourier-transformation the 
spin-orbital interaction Hso readsi^ 



Hso= J2 ["■l,k'-q"-l,-k' 



k.k\q 
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(8) 



where 7 = x,y, z and the V^, ^ are rather lengthy ex- 



X (Vkl-^'A,k+A,-k + Kk'A,k+qC,,k)+h.c 



pressions given in Appendix iDl 

At the points of complete dimerization and in a pure 
spin system the bond-boson description is characterized 
by a vanishing of the inter-dimer coupling. This is not the 
case in the spin-orbital chain, i.e., even at the dimer-point 
Ji, 2=3/4, the matrix elements V^'^, ^ do not vanish and 
the sum of Q and ^ remains a non-trivial interacting 
problem. In turn, spin-triplets can excite orbital-triplets 
by virtue of Hso and vice versa. For an approximate ac- 
count of the resulting renormalization of the one particle 
excitations, we include the effects of Vj!'^, ^ perturbatively 
up to second order in the self energy. The one particle 
excitations result form the poles of the matrix Green's 
function 
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where the second index-pair refers to normal (11 and 22) 
and anomalous (12 and 21) propagators for particles of 
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Figure 2: Self energy of the spin and orbital triplet Green's 
functions, (a) T.a^ y_^{k,uj) (b) (c) Ec,^(fc,t^). Full 

lines denote ^{k,ijj), dashed lines G°^{k,uj) and dotted 



lines G^_, Jfc,w). 



type J^, where v = a^.y^z ov v = Cx,y,z labels spin or or- 
bital triplets. Because G° is calculated from Hlhp it is 
diagonal, i.e. — ^°.2i ~ 0. For the non-diagonal 

elements of Gv this does not apply because the matrix- 
self energy Ei/ is not diagonal. Fig. [2l depicts all contri- 
butions to Sjy up to second order, where each diagram 
occurs with both, an incoming and an outgoing propaga- 
tor on either side, i.e. such that the self energy correc- 
tions can be written in terms of a 2 x 2-matrix notation 
for each of the 3 orbital and 3 spin triplets. These dia- 
grams show that the spin(orbital) triplet motion is renor- 
malized by emission of two orbital(spin) triplets. Earlier 
variational calculations, performed in ref. [13 have been 
based on a subset consisting of exactly these intermediate 
states. We have evaluated the self energy contributions of 
Fig. 121 numerically. The dressed Green's function spec- 
trum displays undamped triplet quasi-particles as well 
as a three-particle continuum. Close to the dimer point 
their respective spectral ranges are well separated. In 
Fig. Elwe show the dispersion of the quasi-particles only. 
For comparison this figure contains complementary data 
form exact diagonalization(ED)i^. Firstly, within a finite 
range of momenta around the zone center and apart from 
a global relative shift of the spectra by ~ 0.1, reasonable 
agreement can be observed for the low-energy excitations, 
i.e. the orbital-triplets. Secondly, with respect to the 
LHP, inclusion of the self-energy from Fig. [^improves the 
agreement with ED. Thirdly, close to the zone boundary 
the ED spectrum refers to the lower edge of the orbital 
two-soliton continuum, rather than the orbital triplet, 
i.e. the two-soliton bound-stateifi. This two-soliton con- 
tinuum is lacking in the approximate bond-operator ap- 
proach. 
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Figure 3: Self energy corrected triplet dispersions of the spin- 
orbital chain. Full lines: spin triplets with (thick line) and 
without (thin line) self-energy corrections. Dashed line: or- 
bital jz-triplets. Dot-dashed line: orbital x,y triplets. The 
dotted line show the border of the multi-soliton continuumi^. 
The symbols denote the spin and orbital triplet energies ob- 
tained by exact diagonalization studies (Lanczos) for a system 
of the length 12 (rungs )^. ED momenta have been back- 
folded into the dimerized Brillouine. 
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Figure 4: Chosen phonon modes, a) transverse distortion 
leading to antiferro-orbital ordering and b) longitudinal dis- 
tortion leading to ferro-orbital ordering. The loops labeled 
orbital one/two denote the 3dyz/3dzx orbital respectively. 
The hatched circles symbolize super-exchange orbitals of p- 
symmetry. Also shown is the static distortion Q. The differ- 
ent size of the orbitals denote the different Coulomb repulsion 
by electrons belonging to neighboring atoms: the larger an or- 
bital is drawn, the lower is its energy. 



III. LATTICE DYNAMICS 

We are now in a position to consider the phonon- 
renormalization due to orbital and spin excitations. To 
this end two steps are required. Firstly, the equilibrium 
structure of the lattice in the orbitally ordered ground 
state has to be clarified. Secondly, the lattice dynam- 
ics has to be evaluated. The first step is specific to 
the material. Here we focus on a situation compatible 
with Na2Ti2Sb20 where the 3dy^/^a;-orbitals form the 
ID chain-structures, shown in Fig. Several static 

Jahn- Teller-type distortions could be envisaged lifting 
the orbital degeneracy. In the remainder of this paper we 
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Figure 5: Self energy of the phonon propagator with (a) 
transverse and (b) longitudinal distortion. As in Fig. [2|full 
lines denote Ga^^ ^{k,Lj), dashed lines Gc^{k,uj) and dotted 
lines Gc^ ^{k,uj). Note, that G"^{k,uj) is calculated with the 

shifted f^^l (see eq. I12|l instead of the pure r'^\ that was 
used in the calculations without any distortion. The explicit 
expressions for the different V^'y ^, i = 1,2 and V^^, j = 1, 4 
are given in appendix IdI 



confine the discussion to the homogeneous ferro{ov 
tudinaiy and antiferro(oT transverse)-0Thit3l static and 
dynamic distortion (c.f. Fig. and b)). The phonon- 
orbiton contribution to the Hamiltonian for both of Fig. 
^) and b) reads 



HpHon = ^ (Q + J3 {bi + bn)) {Tl^ ± Tl^)+u:k,kbibk, ■ 

n 

(10) 

Both phonon modes couple to the z-component of the 
orbital pseudo-spin operator T, where '— ' and '+' refer 
to the transverse, i.e. Fig. ^), and longitudinal, i.e. 
Fig. Eb)) phonon mode. Q accounts for the static Jahn- 
Teller distortion, are phonon destruction(creation) 
operators at site(momentum) n (fc) and ujb,k is the bare 
phonon dispersion. 



operators and after Fourier transformation IjlOII reads: 

k " ' 

+ JsibkTz.k + b-kTz,k + h.c.) +ujb,kblbk , 



TTphon 
trans' 



with a total Hamiltonian oi H ~ H^hp + Hso 
The contribution linear in r^^o resulting from the static 
distortion, i.e. a) in can be eliminated by realizing 

that Hlhp is quadratic in t^^I and by invoking a constant 
shift of the latter operators through: 



^z.k 



QS. 



k.O 



I + J2- J2- 



J2A~ \B) cosk) 



(12) 

On the level of the LHP this does not modify the Tz- 
dispersion. However, inserting (jl2ll into Hso of new 
hopping terms for the spin triplets are generated modi- 
fying their dispersion. Moreover, additional three-triplet 
vertices occur, because of which orbital triplets now cou- 
ple to a two-spin-triplet-continuum. These new terms 
read 



^SO 



^k,q{^l^k^x,y,z,q^x.^y,z,-k-q 

k,q 



^2 —k^^;yj^i~i^^-!y-!^^k-\-q h.c.^ 

+ J2^k,1^^l,kal,y,z,qax.y,z,k+q + h.c.) (13) 
k.q 

They lead to additional self-energy contributions for the 
Green's function of the orbital triplets which have to be 
included in the phononic self-energy because of the cou- 
pling of the orbital triplets to the phonons. Explicit ex- 
pressions for the resulting spin-triplet dispersion (D^^fe on 
the quadratic level, as well as the matrix elements V^!^ 
with i = 1, 2 of this new spin-orbital interaction are given 
in appendix El and appendix ^ respectively. In Fig. El 
we show w^.fc for various values of Q at the dimer point. 
As is obvious, the static Jahn- Teller distortion leads to a 
finite triplet dispersion. 

To evaluate the phonon propagator G, , t {k, u) we use 

the Dyson equation (pj with the self energy depicted in 
Fig. This includes corrections up to second order in 
the phonon-orbiton-spin interactions. The phonon spec- 
trum Sb{k,uj) is then defined as 



Transverse distortion 



First we consider the transverse mode of Fig. ^) , as- 
sociated with antiferro-orbital ordering (AF-00). In this 
case the electrons on two neighboring sites occupy alter- 
nating orbitals. Expressing and r| in terms of bond- 



S'b(fc,u)) = --ImG, f,t(fc,w). 



(14) 



Due to the alternation of the lattice distortion in the 
AF-00 case, i.e. the '+'-sign in there is only a 

linear coupling of the phonon to the orbital z-triplets 
in dJ. Therefore, the dominant renormalization of the 
phonon-spectrum is a direct mixing between z-orbitons 
and phonons. The remaining self-energy corrections 
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Figure 6: Self-energy-corrected spin triplet dispersion without 
(full line) and with different (dashed and dot-dashed lines) 
strength of distortion Q at the dimer point. 
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Figure 7: Phonon spectrum Sb{k,uj) at Ji = 1.2, J2 = 0.8, 
A — B = 1.5, J3 — 0.35 and tot = 1,8 calculated with 
77 = 0.001. The left peak is the one-orbiton excitation peak, 
strongly shifted to lower energies by the distortion of the sys- 
tem. Note the peak position without any distortion (thick 
black line). The position of the phonon-energy peak is shown 
by up-arrows. Because it is an undamped pole of the phonon 
Green's function it is a (5-peak, the height of which is given for 
one example in the figure. Its width is set by the numerical 
broadening rj. The two-particle continuum at uj ~ 2.0 — 2.2 
shows the two van-Hove singularities one expects for a con- 
tinuum of two non-interacting particles. At c^j ~ 2.75 — 3.3 
one finds the three-particle continuum. 
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Figure 8: The phonon peak in the phonon propagator spec- 
trum crossing the orbiton energy at the dimer point, i.e. 
Ji = J2 = 0.75, A ^ B = with = (0.6 + 0.4 cos A:) 
Q — 0.1 and -q = 0.001. The dotted line shows ujb as it is set, 
while the thick full lines denote the peak position of the two 
quasi-particles. The quasi-particle excitations at k ~ n/2 are 
mixtures of orbitons and phonons as expected for interacting 
particles. 



by an arbitrarily small numerical parameter irj. The bare 
phonon energy is momentum independent, i.e. tUb = 1.8. 
This energy remains almost unrenormalized, acquiring a 
very small dispersion and an upward-shift. The energy 
of the orbiton-peak is clearly lowered with respect to its 
bare position (thick solid line in Fig. CJ. This shift in- 
creases with increasing Q, rendering the AF-00 unstable 
at sufHciently large values of the static distortion. Apart 
from the one-particle peaks Fig. [3 displays the antic- 
ipated two- and three-particle continua. Their relative 
weight, as well as the momentum dependence of their 
shape - which exhibits characteristic van-Hove singulari- 
ties - is determined by the matrix elements of ((HJ listed 
in appendix El 

For dispersive phonons, a crossing of the bare phonon 
and orbiton dispersion may occur, where however 
phonon-orbiton coupling will lead to level repulsion in 
the renormaHzed spectra as shown in Fig. |H| In that 
case the quasi-particles in the vicinity of fc = 7r/2 corre- 
spond to an approximately equal-amplitude mixture with 
even and odd parity of orbitons and phonons. 



which are shown in Fig. Ufa) lead to two- and three- 
particle continua. 

Figure H depicts the results of a numerical calcula- 
tion of the phonon-spectrum. In this figure the bare or- 
biton and phonon energies have been chosen to be well 
separated, with a phonon-energy above the orbiton one. 
First, the spectrum displays two one-particle peaks which 
are simple poles of the phonon Green's function, i.e. the 
renormalized phonon and orbiton. Their widths are set 



B. Longitudinal distortion 

Now we consider the longitudinal mode of Fig. E|b), 
which we associate with ferro-orbital ordering (F-00). 
After a Fourier and Bogoliubov transformation and in- 
serting the bond-boson expression for and r| 
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Figure 9: Spectrum of the phonon propagator in the system 
distorted with the longitudinal distortion. Here the parame- 
ters are Ji = 1.2, J2 = 0.8, B = 1.5, J3 = 0.3, = 2.7, 
Q = 0.1 and 77 = 0.0001. Note that there is no feature except 
the phonon-peak at A: = 0. This means this kind of distortion 
will not be seen in spectra of optical experiments. 



reads 

TTphon 

long 



(15) 



k+q 



^y,-k-q'^x,-k 



Again the exact expressions for ^ with i = 3,4 are 
given in appendix El This leads to the diagrams shown 
in Fig. Eb) for the longitudinal distortion. The contribu- 
tion of (I15|l which is static with respect to the distortion 
lifts the degeneracy of the orbital x- and y-triplets. This 
leads to new triplet operators ci{2),k = + (~) 'icx,k + Cy,k 
with an energy-splitting of AlUc = ^Q. Yet, to simplify 
matters the phonon-propagator is evaluated discarding 
this effect, using the Bogoliubov transformed x- and y- 
orbital-triplets of 

While for the AF-00 a linear coupling of phonons 
and orbitons results from lfT?Hl . the '—'-sign in the lat- 
ter eqn. leads to ifTfill . i.e. the decay of phonons into 
orbital triplet-pairs in the case of F-00. In turn the pri- 
mary effect on the phonon spectrum in the F-00 is the 
appearance of the a two-orbiton continuum. This can be 
observed in Fig. |2| which depicts the numerical solution 
of the Dyson equation of Fig. |3b). The spectrum shows 
a renormalized simple phonon-pole the width of which 
is set by a numerical broadening of irj = ilO~^. The 
weak phonon-dispersion observable is an interaction ef- 
fects, since the bare phonon has been chosen to be purely 
optical. In addition to the renormalized phonon Fig. |2| 
exhibits a two-orbiton continuum, the spectral weight of 
which is more than one order of magnitude larger than for 
the AF-00 case. In the case of the finite splitting Awc 



included Fig. Elwould show three superimposed continua, 
a case which we leave to future analysis. As for the AF- 
00 the continuum shows characteristic van-Hove type 
singularities. Most important however, at momentum 
k = the spectrum resembles a pure phonon excitation 
only. This is caused by a vanishing of the matrix element 
^k.q at g = (see App. El and impHes that this kind of 
distortion will not have any effect on spectra of optical 
experiments e.g. Raman-scattering. 



IV. CONCLUSION 

In conclusion we have studied the interplay of spin, or- 
bit and lattice dynamics in a model of an anisotropic spin- 
orbital coupled chain in the vicinity of the dimer-point. 
The spin-orbital sector of this model has been treated 
by mapping onto interacting bond-bosons with a singlet 
condensate and massive triplet modes. We have analyzed 
the lifting of the orbital degeneracy for two kinds of dy- 
namic Jahn- Teller distortions leading either to antiferro- 
or ferro-orbital ordering. Studying the phonon dynam- 
ics in both cases we have found shake-up sidebands which 
are critically dependent on the type of distortion. In par- 
ticular we have found (no) optical activity of the phonon 
in the transverse (longitudinal) distortion. 
Acknowledgment: We thank C. Jurecka for stimulat- 
ing discussions in the beginning of this work. This re- 
search was supported in part by the DFG, Schwerpunk- 
tprogramm 1073. 



Appendix A: MATRIX ELEMENTS OF Ho 

For the matrix elements of the quadratic part of the 
Hamiltonian (Eq. ^ one finds M,j^„ = Ji, = 
,h{l + A/2), Af56„ = J2, Af,i^i,„ = Af2^i,„ = Mi^;^„ = 
-l/4(Ji - 3/4 - 1/4B), „ = Af,f_ i,„ = 

Mt\n = -1/4(^2 - 3/4) and A/^e ^ „ = 
Mf„i,„ = Aff_i,„ - Ml\„ - -l/4(J2(l + A)-3/4(l + 
B)). All the other matrix elements are zero. 



Appendix B: SPIN AND ORBITAL 
TRIPLET-DISPERSIONS AND BOGOLIUBOV 
COEFFICIENTS 

Here the spin and orbital triplet dispersions after LHP 
approximation (0 are given. 



UJa,k 



Jl 



1- ( 1-^(1 + 5) )cosfc 
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1 I cos k 
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1 (l + B)+A]cosk 
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(Bl) 



7 



The Bogoliubov coefficients are defined as 



2 

9k 



cos k 



•h - i(J2(l + A) - 1(1 + B))cosk 



J2(l + 4)-^(J2-|)cOSfc 



(B2) 



and vl^ul- 1, = .g2 - 1, = ^2 _ 



Appendix C: SPIN-TRIPLET DISPERSION FOR 
STATIC TRANSVERSAL DISTORTION 



For tfie spin-triplet dispersion Cj^^k for tlie static 
transversal distortion one finds: 



+ {^^,k 



64g2(4(l + B)ukVk + [uk + vkf cos fc) 
(7 + 35 - 4AJ2)4 

32(2(1 + B)(u| + vl) + [uk + Vkf cos fc) 



(7 + 35 - 4AJ2)2 



(CI) 



where u)-y^k is the spin triplet dispersion of the undistorted 
model and and Vk are the Bogoliubov-coefficients (Ap- 
pendix El of the transformations in 



^k,k',q 



y3,z 
^k,k\q 



^kM',q 



= -^[Tk'qIk+q,kWk'-qUJk' 
O 

+ '^BiT^_gUk'-qVk'ik+q-ik] (D3) 

= --^[B2Tk'qOk+q,kWk'-qWk' 

+ 2T^_^Uk'-qVk'0k+q0k] (D4) 

= --^[2Tk-qgk'-qhk'Wk+qWk 

+ BiT*,^^Wk+q,kik'-qik'] (D5) 

= -\[B22Tk-qmk'^qnk'Wk+qWk 

o 

+ T* ,^M^fc+,,feOfc._gOfc.] (D6) 

+ Tk',q{uk'-qVk' + Vk'-qUk')Ik+q.k 

+ B,Tl_^Wk'-q,k'{gk+qhk + hk+qQk)] (D7) 

= -\[{Tk',qB2 + Tl_^)Ok + q,kWk',k'~q 

+ B2Tk'^g{Uk'^qVk' +Vk'-qUk')Ok+q,k 

+ T^_qWk'-q.k'{rnk+qnk+nk+qmk)]. (D8) 



The matrix elements of the spin-orbital interaction Hso 
arising from the shifted orbital triplets are: 



Appendix D: MATRIX ELEMENTS OF THE 
SPIN-ORBITAL INTERACTION Hso 



For completeness the matrix elements of the spin or- 
bital interaction ^ are listed. Here T^^y = e~^^ +e'^^^~y\ 
With Uk + Vk = Wk, gk + hk = ik, ruk + Uk = Ok, 
UkUk' + VkVk' = Wk,k', gkgk' + hkhk' = Ik,k', mkruk' + 
rikUk' = Ok,k' and 1 -h B = 5i, 1 B/2 = B2: 



^k,k',q 



^k,k',q 



-^[{2Tk'^q + 2BlT^_g)gk+qhkUk'-qVk' 

Tk',q{gk+qhkUk'-qUk' + hk+qgkVk' -qVk') 
BlT^_q[gk+qgkUk' -qVk' + 

Uk'Vk'-qhk+qhk)] (Dl) 
-^[{2B2Tk'^q + 2T^_g)mk+qnkUk'-qVk' 

B2Tk' ,q{mk+qnkUk' -qUk' + nk+qirikVk'-qVk') 
Tk^_g{mk+qmkUk' -qVk' + 

Uk'Vk'-qUk+qnk)] (D2) 



-FY,[<^osq- fk-{Wk^k+q)+ 



k,q 



k,q 



k.q 

+ (cosg + cos(fc + q) + 2(1 + B) ■ [l + cosfc)) x 

X {gkUqVk+q + hkVqUk+q)] (D9) 

^ [2 cos q- {ik) ■ {ukVk+q + VqUk+q) + 
k.q 

+ (cos q + cos(fc 4- g) • 2 • (1 + B) • (1 + COS fc)) X 

X (ik) ■ {Wk,k + q)] (DIO) 

(Dll) 
(D12) 



'2tOk. 



k+q 



Vk^q = 2i{mknk+q - rrik+qUk) 



where F = Q/{{J2 + 4X)uJbfi + 4J|), X ^ -1/4(J2 - 
3/4 + J2A - 3/4i?) and Uk, Vk, gk, hk, ruk, Uk are the 
Bogoliubov-coefficients (Appendix El of the transforma- 
tions in ©. 
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